Spin polarized transport in the weak link between f-wave superconductors 
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The spin current in the Josephson junction as a weak-link (interface) between misorientated triplet 
superconductors is investigated theoretically for the models of the order parameter in UPts. Green 
functions of the system are obtained from the quasiclassiacl Eilenberger equations. The analytical 
results for the charge and spin currents are illustrated by numerical calculations for the certain 
misorientation angles of gap vector of superconductors. As the main result of this paper, it is found 
that, at some values of the phase difference, at which the charge current is exactly zero, the spin 
current has its maximum value. Furthermore, it is shown that the origin of spin current is the 
misorientation between gap vectors of triplet superconductors. 

PACS numbers: 74.50.+r, 74.70.Pq, 74.70.Tx, 72.25.-b 



INTRODUCTION 



Triplet superconductivity has become one of the most 
interesting topics of condensed matter physics 1 ' 2 , partic- 
ularly in view of the recently discovered ferromagnetic 
superconductivity 3 -^. The mechanism of pairing, physics 
of interaction and gap structure in this type of super- 
conductors have been the subject of many experimen- 
tal and theoretical works^^ ^ 1 ' 1 1 . The Cooper pair- 
ing in the triplet superconductors has been investigated, 
for example, using the thermal conductivity in Refsj 12 ' 13 
and Knight shift experiments in Refs i 14 ' 15 . Also, the 
Josephson effect in the point contact between triplet su- 
perconductors has been studied in Ref. 16 . These weak- 
link structures have been used to demonstrate the or- 
der parameter symmetry in RefJi. Eventually, the /- 
wave symmetry of order parameter has been proposed 
for UPt^ and SV2-RUO4 compounds. In addition, the 
spin polarized transport through the systems consist- 
ing of superconductors, normal metals, ferromagnetic 
layers and other structures as one of the modern top- 
ics of mesoscopic physics, has attracted much attention 
recentl y 18 ' 19 i 20 i 21 . In the present paper, the ballistic 
Josephson weak-link as the interface between two bulk of 
/-wave superconductors with different orientations of the 
crystallographic axes has been investigated. It is shown 
that the current-phase dependencies are totally differ- 
ent from the current-phase dependencies of the junction 
between conventional (s-wave) superconductor s 22 ' 23 and 
high T c (d-wave) superconductors 2 ^. It is found that for 
the certain values of the misorientation, the spin-current 
in the both directions, tangential and perpendicular to 
the interface, may exist and it has totally unusual de- 
pendence on the external phase difference. The effect of 
misorientation on the spin current is investigated. It is 
observed that the misorientation between gap vectors is 
the origin of the spin current. As the important result of 
this paper, it is obtained that, at some of certain values 
of phase difference, at which the charge current is zero, 
the spin current has the finite value. Another result of 



the paper is the capability of this proposed experiment 
for polarization of the spin transport using the junction 
between /-wave superconductors. Eventually, one of the 
states and geometries of our system can be used as a 
switch which is able to divide the spin and charge cur- 
rents into two parts: parallel and perpendicular to the 
interface. 

The organization of the rest of this paper is as follows. 
In Sec. II we describe our configuration, which is investi- 
gated. For a non-self-consistent model of the order pa- 
rameter, the quasiclassiacl Eilenberger equations 3 -! are 
solved and suitable Green functions are obtained analyt- 
ically. In Sec. Ill the obtained formulas for the Green 
functions are used for calculation the charge and spin cur- 
rent densities at the interface. An analysis of numerical 
results will be done in Sec. IV. The paper will be finished 
with some conclusions in Sec.V. 



II. FORMALISM AND BASIC EQUATIONS 

We consider a model of a flat interface y = between 
two misorientated /—wave superconducting half-spaces 
(Fig.l) as a ballistic Josephson junction. In the quasiclas- 
sical ballistic approach, in order to calculate the charge 
and spin current, we use "transport-like" equations 3 ^ for 
the energy integrated Green functions g(yp,T, e m ) 



s m a 3 + iA,g =0 



and the normalization condition 
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where e m = 
m = 0,1,2. 



7rT(2m + 1) are discrete Matsubara energies 
.., T is the temperature and Vf is the Fermi 



velocity and 03 — 03® I in which aj (j — 1, 2, 3) are Pauli 
matrices. The Matsubara propagator g can be written in 
the standard form: 



91 +glCT 

i&2 (53 + g3<5-) 



(52 

94 - 



-g 2 erj 102 

<7 2 g40'd'2 



(3) 



2 




where, the matrix structure of the off-diagonal self en- 
ergy A in the Nambu space is 

In this paper, the unitary states, for which d x d* = 0, 
is investigated. Also, the unitary states vectors di2 can 
be written as 

d n = A n exp iip n , (5) 

where A1.2 are the real vectors in the left and right sides 
of the junction. The gap (order parameter) vector d has 
to be determined from the self-consistency equation: 

d (v F , r) = irTN (0) £ (V (v F , v F ) g 2 (v' F , r, e m )) (6) 

m 

where V (vf,v' f ), is a potential of pairing interaction, 
(...) stands for averaging over the directions of an elec- 
tron momentum on the Fermi surface v' F and N (0) is 
the electron density of states at the Fermi level of en- 
ergy. Solutions to Eqs. (fTJ) and © must satisfy the 
conditions for Green functions and vector d in the bulks 
of the superconductors far from the interface as follow: 

„ , , x s m a 3 + iA 21 

9 (±00) = / =; (7) 

\M + |d 2 ,i| 2 
d(±oo) = d 2 , 1 (v F )exp^Ty) , (8) 

where <fi is the external phase difference between the order 
parameters of the bulks. Eqs. flU and © have to be sup- 
plemented by the continuity conditions at the interface 
between superconductors. For all quasiparticle trajecto- 
ries, the Green functions satisfy the boundary conditions 
both in the right and left bulks as well as at the interface. 



FIG. 2: Charge and spin current (s z ) versus the phase dif- 
ference <j> for the planar state (|16[) . geometry (i) and the 
different misorientations. Currents are given in units of 
jo,e = f eN{0)v F A (0) and j 0>s = f hN(0)v F A (0) respec- 
tively. 

The system of equations (JTJand @ can be solved only nu- 
merically. For unconventional superconductors such solu- 
tion requires the information of the function V (vjr, v F ). 
This information, as that of the nature of unconven- 
tional superconductivity in novel compounds, in most 
cases is unknown. Usually, the spatial variation of the 
order parameter and its dependence on the momentum 
direction can be separated in the form of A(v>,y) = 
A(vf)\I>(j/). It has been shown that the absolute value 
of a self-consistent order parameter and ^(y) are sup- 
pressed near the interface and at the distances of the 
order of the coherence length, while its dependence on 
the direction in the momentum space (A(vp)) remains 
unaltered™. Consequently, this suppression doesn't in- 
fluence the Josephson effect drastically. This suppression 
of the order parameter keeps the current-phase depen- 
dence unchanged but, it changes the amplitude value of 
the current. For example, it has been verified in Refs^ 
for the junction between unconventional d-wave, in Ref. 25 
for the case of "/-wave" superconductors and in Rcf. 26 
for pinholes in 3 He that, there is a good qualitative agree- 
ment between self-consistent and non-self-consistent re- 
sults. Also, it has been observed that the results of the 
non-self-consistent investigation oiD — N—D structure in 
the paper— are coincident with the experimental results 
of the paper— and the results of the non-self-consistent 
model in paper-2^ are similar to the superfluid weak-link 
experiment^,. Note that any self-consistent solution can 
be found for some models of pairing potential and the 
certain values of parameters. Consequently, despite the 
fact that this estimation cannot be applied directly for 
a quantitative analyze of the real experiment, only a 
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FIG. 3: Charge and spin current (s y ) versus the phase differ- 
ence 4> for the axial state (|15p . geometry (ii) and the different 
misorientations (y-component). 



FIG. 4: Charge and spin current (s y ) versus the phase differ- 
ence 4> f° r the axial state (|15p . geometry (ii) and the different 
misorientations (x-component) 



qualitative comparison of calculated and experimental 
current-phase relations is possible. In our calculations, 
a simple model of the constant order parameter up to 
the interface is considered and the pair breaking and the 
scattering on the interface are ignored. We believe that 
under these strong assumptions our results describe the 
real situation qualitatively. In the framework of such 
model, the analytical expressions for the charge and spin 
current can be obtained for an arbitrary form of the order 
parameter. 



III. ANALYTICAL RESULTS. 

The solution of Eqs. ((T|) and ^ allows us to calculate 
the charge and spin current densities. The expression for 
the charge current is: 



Eilenberger equations |1} for Green functions g, which 
are supplemented by the condition of continuity of solu- 
tions across the interface, y = 0, and the boundary con- 
ditions at the bulks, are solved for a non-self-consistent 
model of the order parameter analytically. Two diago- 
nal terms of Green matrix which determine the current 
densities at the interface, y = 0, are following. For the 
relative term to the charge current we obtain: 



9i (0) 



£ m (»i +Sl 2 )cos^ + i?7sin/3(Sli»2 +Jm) 
ir)smj3e m (^i + ^2) + cos/3(f7i^ 2 + e 2 rn ) + AiA 2 

and for the case of spin current we have: 



Si (0) 



11 Al X A 2 2 7 [(B-l) 2 cxp(z/3)(7 ? f7 1 +e m )x 



j e (r) = 2ineTN (0) J2 ( v ^5i (vf, r, e m )) (9) (r^ 2 + e m ) — (B + l) 2 exp(-^)(r^ 2 - e m ){r)Q,i - e m )| 



and for the spin current we have: 



j Si (r) = 2 l ^(-)T7V(0)^(v F (e lgl (v F ,r, £ „ l ))) (10) B = 



where, e^= (x, y,z). We assume that the order parame- 
ter does not depend on coordinates and in each half-space 
it equals to its value ([5]) far form the interface in the 
left or right bulks. For such a model, the current-phase 
dependence of a Josephson junction can be calculated 
analytically. It enables us to analyze the main features 
of current-phase dependence for the different models of 
the order parameter of "/-wave" superconductivity. The 



(12) 

where 77 = sgn{v y ), fl„ = \je 2 m + |d„| 2 , f3 = ipi-ip2+<f>, 

??e m (£li + tt 2 ) cos (3 + zsin/?(»i» 2 + e 2 J 

M7sin/3£ m (fii + fi 2 ) +cos/3(0 1 fJ 2 + £ 2 J + A i A 2 

(13) 



and 



A = 



A 1 A 2 (B- l)exp(z/3) + A^jB + 1) cxp(-i/3) 

(77^1 - £ m )(T)Q2 ~ Em) + £ m ){V^2 + £m) 

(14) 

Also, n = 1, 2 label the left and right half-spaces respec- 
tively. We consider a rotation R only in the right super- 
conductor (see, Fig.l), (i.e., d 2 (k) = J^di(JS _1 k), k is 
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FIG. 5: Charge and spin current (s y ) versus the phase differ- 
ence (j> for the axial state (|15p . geometry (ii) and the different 
misorientations (z-component) 

the unit vector in the momentum space). The crystallo- 
graphic c-axis in the left half-space is selected parallel to 
the partition between the superconductors (along z-axis 
in Fig.l). To illustrate the results obtained by computing 
the formula (111|12|) . we plot the current-phase diagrams 
for the different models of the "/-wave" pairing symme- 
try (| 1 5I16[) and for two different geometries. These two 
geometries are corresponding to the different orientations 
of the crystals in the right and left sides of the interface 
(see, Fig.l): 

(i) The basal a6-plane in the right side has been rotated 
around the c-axis by a; Ci||c2. 

(ii) The c-axis in the right side has been rotated around 
the axis perpendicular to the interface (y-axis in Fig.l) 
by a; bi||b 2 . 

Further calculations require a certain model of the gap 
vector (vector of order parameter) d. 

IV. ANALYSIS OF NUMERICAL RESULTS 

In the present paper, two most probable forms of the /- 
wave order parameter vector in UPt^ are considered. The 
first model which is successful to explain the properties 
of the B-phase of UPt^ is the axial state. This sate 
describes the strong spin-orbital coupling with vector d 
directed along the c axis of the lattice and it is: 

d(T,v F ) = A (T)zk z (k x + iky f . (15) 

The coordinate axes x, y, z here and below are chosen 
along the crystallographic axes a, b, c in the left side of 
Figffl The function A = A (T) in Eq.|T5D and below 
describes the dependence of the order parameter d on 



FIG. 6: Tangential spin current (s z ) versus the phase differ- 
ence 4> for the planar state (|16[) . geometry (ii) and the different 
misorientations. The perpendicular component (y- direct ion) 
of the spin current is absent. 

the temperature T (our numerical calculations have been 
done at the temperatures close to the T — 0). The second 
model of the order parameter which describes the weak 
spin-orbital coupling in UPt% states, is the unitary planar 
state. The planar model of gap vector is: 

d(T, v F ) = A (T)fc z (x (k 2 x - k 2 y ) + y2k x k y ). (16) 

Using these two models of order parameters (|15I16|) and 
solutions to the Eilenberger equations (|llll2p , we have 
calculated the spin current and charge current densities 
at the interface numerically. These numerical results are 
listed below: 

1) The spin current can be present, only when misorien- 
tation between gap vectors exists. Because in our Green 
function (|12[) . the spin current is proportional to the 
"cross product" between the left and right gap vectors. 
For instance, the spin current for the case of the axial 
state (|TS1) and geometry (i) is zero, because both of the 
gap vectors are in the same direction (z). (Geometry (i) 
is a rotation as much as a, around the z axis). 

2) In Fig[2]it is shown that for the planar state and geom- 
etry (i), it is possible to observe the current of s z in the 
direction perpendicular to the interface, but in Figs l3l4l 
and it is demonstrated that, for the axial state and 
geometry (ii), only the current of s y can be observed. 
Consequently, this kind of junction can be applied as a 
polarizer or filter for the spin currents. 

However, for the planar state and geometry (ii), all terms 
of the spin current (s x , s v and s z ) can be observed (see 
EqmD. 

3) In Figs[3][5] (planar states), it is shown that the value 
of the phase differences in which the currents are in the 
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FIG. 7: The perpendicular component of the charge current 
(Josephson) versus the phase difference <j> for the planar state 
(|16|) . geometry (ii) and the different misorientations. The 
tangential components (x and 2-directions) are absent. 



FIG. 8: Tangential spin currents (s x ) versus the phase differ- 
ence 4> for the planar state (JT6J) , geometry (ii) and the different 
misorientations. The perpendicular component (y- direct ion) 
of the spin current is absent. 



maxima, minima and zero values, are not very sensitive to 
the misorientation angle a, while the amplitude of max- 
ima and minima, are strongly dependent on the value of 
misorientation a. 

4) In the Figs l2l3l and Fig[5J[71 while the charge currents 
are the odd functions of <f> with respect to the line of 
<f> — iv, the spin currents are even functions of the phase 
difference; j e (</> = n + 8<j>) — — j e (</> — tt ~ 8(f) and for the 
spin current ] Si {4> = 7r + 8<f>) = j Si (4> = tt — 8(f)). On the 
contrary, in the Figsf4]and[5l the charge and spin currents 
are even and odd functions of 4> with respect to the line 
of (f) — 7T, respectively; j e ((j> = ir + 8(f) = j e {4> = 71 ~ ^4>) 
and j Si (4> = 7r + 84>) = —j ai (<f> = tt - 8(f). 

5) In Fig[21 the perpendicular component of the spin and 
charge current in terms of the external phase difference <p 
for the case of the planar state geometry (i) and for 
two different misorientations are plotted. The solid line 
is the charge current-phase dependence^. Also, at the 
<j> = 0, (f> = 7r and (j) = 27r, the charge current (Joseph- 
son current) is zero while the spin current has the finite 
value. 

6) The perpendicular component of the charge (Joseph- 
son current) and spin current for the case of the axial 
state (fL3|) and geometry (ii) are plotted in Fig[31 and the 
tangential components of them, are plotted in Figs l4l5l 
The charge current-phase diagrams have been obtained 
before in paper— and they are totally different from the 
case of conventional superconductors in the paper—. At 
the phase values of <p = 0, <p = tt and (j> = 2ir, in which 
the charge current is exactly zero, the spin current has 
the finite values and may select its maximum value. In 
Figs l2l3l and specially FigJSl for a small value of misori- 



entation we have a very long but narrow peak in the spin 
current phase diagram, close to the tf> = 7r. 
7) Both the planar state with geometry (i) and the axial 
state with geometry (ii) can be applied as the filter for 
polarization of the spin transport (see Fig^2][5]), the for- 
mer transports only the s z but the latter case flows the 
s y (see rT2"T) . In addition, the planar states with geome- 
try (ii) can be used as a switch for the spin and charge 
current into two directions: parallel and perpendicular to 
the interface. In this case, the spin and charge currents 
select only one of the directions parallel or perpendicu- 
lar to the interface. Namely, it is impossible to observe 
the tangential and perpendicular components of the cur- 
rents at the same time for planar state with geometry 
(ii)(Fi g sEU). 



V. CONCLUSIONS 

We have theoretically studied the spin current in the 
ballistic Josephson junction in the model of an ideal 
transparent interface between two misorientated /-wave 
superconductors which are subject to a phase difference 
<p. Our analysis has shown that the misorientation and 
different models of the gap vectors influence the spin cur- 
rent. This has been shown for the charge current in the 
paper—. The misorientation changes strongly the critical 
values of both the spin current and charge current. It has 
been obtained that the spin current is the result of the 
misorientation between the gap vectors. Furthermore, it 
is observed that the different models of the gap vectors 
and geometries can be applied to the polarization of the 
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FIG. 9: Perpendicular component of the spin current (s y ) 
versus the phase difference <f) for the planar state (|16|) . ge- 
ometry (ii) and the different misorientations. The tangential 
components (x and z directions) are absent. 

spin transport. Another result of these calculations is the 
state in which the currents select one of the two possible 
directions ( perpendicular and parallel to the interface) 



to flow. This property can be used as a switch to con- 
trol the direction of the charge and spin current. Finally, 
as an interesting and new result, it is observed that at 
some certain values of the phase difference 0, the charge- 
current vanishes while the spin-current flows, although 
the carriers of both spin and charge are the same (elec- 
trons). The spatial variation of the phase of the order 
parameter plays a role as the origin of the charge current 
and, similarly, due to the broken Q s P m ~ orb it symmetry, 
a spatial difference of the gap vectors in two half-spaces 
causes spin currents. This is because there is a position- 
dependent phase difference between "spin up" and "spin 
down" Cooper pairs and, although the total charge cur- 
rent vanishes, there can be a net transfer of the spin. 
Therefore, in our system, there is a discontinuous jump 
between the gap vectors and, consequently the spin cur- 
rents should generally be present. For instance, if up-spin 
states and down-spin states have a velocity in the op- 
posite direction, the charge currents cancel each other 
whereas the spin current is being transported. Math- 
ematically speaking, j c har ge = h + ji,js P in = Jt - Jl, so 
it is possible to find the state in which one of these current 
terms is zero and the other term has the finite valued, 
fn addition the spin imbalance which is the result of the 
different density of states for "spin-up" and "spin-down" 
can be the other reason of spin current^. In conclusion, 
the spin current in the absence of the charge current can 
be observed. 
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